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Abstract 

A generalized gauge invariant Ising model on random surfaces with 
non-trivial topology is proposed and investigated with the dual trans- 
formation. It is proved that the model is self-dual in case of a self-dual 
lattice. In special cases the model reduces to the known solvable Ising- 
type models. 
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1 Introduction 



In recent years great interest has turned to the non-critical string theory with 
central charge C < 1, which was claimed to be non-perturbatively solvable 
by relating its dual model to certain matrix representations. 

The string theory with C = 1/2 can be viewed as the Ising model on ran- 
dom surfaces. In the non-perturbative calculation, the continuum surface is 
tessellated by polygons with p edges. Then through the dual transformation, 
the Ising model is related to its dual (Ising) model on the 0^ lattice. The 
latter has a two-matrix representation which is non-perturbatively solvable 
around the double scaling point 0, |^ . 

The dual tranformation, however, is not straightforward for non-planar 
lattices because the high-temperature expansion contains non-boundary one- 
cycles. As early as 1971, Wegner ^ pointed out that the traditional duality 
relation is valid only if the completeness condition is fulfilled. For two di- 
mensional lattices Wegner's completeness condition means nothing but the 
planar lattice. 

Recently a generalized duality relation valid for both planar and non- 
planar lattices has been obtained The remarkable fact is that the 
dual model of an Ising model is not simply an Ising model, as one implicitly 
supposed in the two- matrix representation of the string theory with C = 1/2. 
It has also been shown that the model, whose dual model is an Ising model 
coupled to an external field, turns out to be a gauge invariant Ising model 
where non-boundary one-cycles do not appear due to the gauge invariance. 
This was first suggested by Wegner and Balian et al. for a regular 
planar lattice. For what has recently been done, the exact statement should, 
therefore, be that it is the gauge invariant Ising model on random surfaces 
that has the solvable two-matrix representation. Whether the gauge invariant 
Ising model and the Ising model coupled to an external field belong to the 
same universality class remains open.| 

All these indicate that the gauge field in two dimensions is far from trivial, 
especially in the case of non-planar lattices. It is therefore interesting to 
investigate Ising-type models on random surfaces with gauge interaction and 
to study the topological effects of the gauge field. 

^ One may directly discretize the continuum theory by putting the spins on the dual 
lattice, but in this case one also has to know whether this discretized model is in the same 
universality class of the original one or not. 
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On the other hand, as is well known, the self-duality is one of the most 
attractive points of the Ising model, which help us to understand the phase 
structure. The self-duality is in general valid only when the lattice is self-dual. 
For the planar random surfaces one may get the self-duality by summing up 
all the lattices. This has recently been discussed 0. On non-planar lattices 
the situation is more complicated. The self-duality is spoiled by the non- 
trivial topology of the lattice even if the lattice is self-dual. What is then a 
self- dual model on non-planar random surfaces? 

Furthermore the dual transformation itself is also useful in studying the 
relation between two Ising-type models |p. In two dimensions the topo- 
logical properties of the surfaces are closely related to boundary conditions. 
Therefore it is important to know how the dual transformation goes for topo- 
logically different lattices. 

In order to tackle these problems, in this paper the simple gauge in- 
variant Ising model |^ is generalized such that two Ising-spins interact 
with a gauge field. It will be seen that this model is self-dual even for a 
non-planar lattice if the lattice is self-dual, because the gauge field kills the 
non-boundary one-cycles in the high-temperature representation. This is the 
simplest self-dual Ising-type model on the lattice with non-trivial topology. 
Under certain limits the model reduces to relatively simpler ones, which are 
already known to be physically interesting and important. We can demon- 
strate the calculation of Wilson loops and show how the topology plays its 
role. 

In the next section the dual transformation will be carried out. Several 
limiting cases are considered in section 3. Finally some discussions follow. 

2 A self-dual model 

Let us consider a system on a random lattice L with two sets of Ising spins 
and a Z2 gauge field, whose partition function is 

z = \ Y [JUijSiSj + j'UijS-s'j^ 

{s.,s'^,u,j} [<ij>eL^ 

+ Y PUm+Y ^S,Sl 
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where L°,L^ and are respectively the set of sites, bonds and plaquettes; 
Si = ±1 and S^ = ±1 are the spin variables on a site, Uij = ±1 is the gauge 
field on a bond and Um is the plaquette variable defined as the product of 
the gauge fields along the boundary of the plaquette. This model looks like 
the 'spin-Schwinger model'. J^ieLO XSiS^ corresponds to the 'mass' term. 

Here we should mention, as in the case of the simple gauge invariant Ising 
model [0], in the 'vacuum sector' the spins {S^} can be gauged away by the 
following tranformation 

Uij = UijSj^Sj, Si = SiS^ (2) 

. However, for the convenience of discussions we will keep it because it 
preserves explicitly the gauge invariance of the system. 

For a given random lattice L, one can construct a dual lattice L^, and a 
model on it with the partition function 



(3) 



Zd = E *^^P ) E [jDVmnTmTn + JD^m 

+ J2 l^oVi + E ^DTmT^ 

where L'j^jLjj and Lj^ are respectively the set of dual sites, dual bonds and 
dual plaquettes; r^, t„ = ±1 are dual spin variables on dual sites, Vmn = ±1 
is the dual gauge field on a dual bond and Vi is the dual plaquette variable 
defined as the product of the dual gauge fields along the boundary of the 
dual plaquette. 

Following the standard procedure for the dual transformation, we will 
prove that the model with Z is dual to the one with Z^. Rigorously speaking, 

Z = cZn, c = 23(^-^)/2(sinh2Jsinh 2j')^/2(gi^j^2/3)^/2(gjnl^2A)^/2 (4) 

if 

= -\n (coth A) , Ad = - In (coth p) (5) 



and either 



= - In (coth J), = - In (coth J ) (6) 
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Jd = - In (coth J ), = - In (coth J), (7) 

where B, N and P are respectively the numbers of the bonds, sites and 
plaqettes on the original lattice. 

Note that if L Ld, Z and Z^ are respectively partition functions for 
two different models, which are dual each other. But in case of L = L^, these 
two model become the same, i.e. we have a self-dual model. Then the dual 
relations in (^0) may give some insight in the phase structure. From the 
dual relation @ we know that the critical points appear always by pair in 
the (/?, A) plane. It is, however, not possible to locate the critical point even 
when it is assumed to be unique, as in the case of the simple Ising model. 
From the dual relations (§[-0) one can see that the critical points are grouped 
in four in the {J, J ) plane. If there is only a unique critical point in the 
{J, J ) plane, it is {J*, J*) with J* = | In (coth J*). To make full use of the 
dual relations further investigations are needed. 

In order to get the duality relations (§-0), we simply write down the high- 
temperature representation for Z and the low-temperature representation for 
Zd and compare the two expressions. 

2.1 The high-temperature representation of Z 

Noting that Si, S*- and Uij all take values ±1, the partition function Z in 
Eq.(|I|), can be written as 



Z = ^ Y[ (cosh J + UijSiSj sinh J) Y]_ (cosh J + UijS^S^ sinh J ^ 

Y[ (cosh A + SiSi sinh A^ exp < ^ (3Um \ 

(8) 

The product Y[<ij>eL^{cosh.J + UySiSj sinh J) can be expanded into 2^ 
terms. For a given term, a bond < ij > contributes a factor of either 
cosh J or UijSiSj sinh J. Let us mark the bond in case of UijSiSj sinh J. 
All the marked bonds for this given term assemble as a one-chain 7 on the 
lattice. It is easy to see that this is a one-to-one correspondence between all 
the terms in the expansion and all the one-chains on the lattice. Similarly 
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7 + 7' = 70 



Figure 1: This is a handle of the random surface. The 70 on the right side 
is a boundary one-cycle. The 70 on the left is a non-boundary one-cycle 

n<jj>gLi (cosh j' + UijS[S'j sinh j') can be expanded according to another set 
of all one-chains 7'. For niGLo(coshA + ^jSj- sinh A) we mark the site i if 
a term in the expansion picks up SiS^ sinh A. All the marked sites in the 
term form a zero-chain 9. This is also a one-to-one correspondence. After 
performing the summation over {Si} and {5*^}, we get 

Z = 2'''Y. n sinhJ n coshJ \{ sinh/ 

{C^ijlT^'eCi eeco <ii>e7 <i3><^'i <ij>e-y' 

Y[ cosh J sinh A J]^ cosh A Y\_ ^ij 6xp < (3U, 

where and are the sets of zero-chains and one-chains respectively, 

7o = 7 + 7', and 

e + d-f = 0, 9 + d-f' =0 (10) 

Here d'y and d'y' denote respectively the boundaries of 7 and7'. The equation 
( p!0| ) simply means that 7 and 7 are connected by 9 resulting that 70 must 
be a one-cycle. 

Here we should stress that up to now the gauge invariance has not been 
considered and 70 can be either a boundary one-cycle, which circles a certain 
area of the surface, or a non-boundary one-cycle, which bounds no area. This 
is shown in Fig |I|. 

Similarly one can carry out the summation over gauge fields. If 70 is a 
boundary one-cycle. 



{Uij} <ij>ejo [mGL2 

2^1 n sinh/3 cosh/5+ [] cosh/3 sinh/3 

[^mGS-yo mfS'fo mG<57o m^S'yo 



(9) 
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Figure 2: This is a part of the original lattice and the dual lattice. Here 
the circles and the dotted lines represent the dual sites and dual bonds. The 
real dots and lines are the sites and bonds on the original lattice. Vmn is the 
dual gauge field. 

where ^70 is the set of all plaquettes with boundary 70. The two terms on the 
right-hand side represent the two ways to fill up 70, 'inside' and 'outside'. If 
7o is a non-boundary one-cycle, the summation over gauge fields gives zero. 
Therefore 



Z = 2^+2iv^ ^ j-j ^.^^j Yl cosh J 

Y[ sinh J Y\_ cosh J Y]_ sinh A cosh A 

<jj>e7+70 <«j>^7+70 «e97 i^d-y (12) 

Yl sinh/? Yl cosh f3+ J| cosh/3 JJ sinh/3 

where Fq is the set of all boundary one-cycles. This is the so-called high- 
temperature representation of Z. The gauge field damps the non- boundary 
one-cycles induced by the non-trivial topology of the lattice. 

2.2 The low-temperature representation 

The configurations of the dual gauge fields can graphically be represented by 
all the one-chains on the original lattice. Let us denote the dual bond of 
< mn > by < mn >*. For a given configuration, the corresponding one-chain 
7 is obtained as follows: if Vmn = —1, draw a line on the bond < mn >* 
on the original lattice; if Vmn = 1, do nothing. It is not difficult to see that 
on the boundary of 7 the dual plaquette variables are —1 and otherwise 
This is shown in Fig ^ 
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Figure 3: (a) and (b) are the two configurations of dual spins {t^} corre- 
sponding to tlie boundary one-cycle 70. 

Therefore can be written as 

^D= n exp {-JoT^Tn - JDr'rrJr] 

{r™,T;„}7eCi <mn>*G7 

W exp I JDTmr„ + JdvI}^ 

<mn>* ^7 

n e-^- n E adt^t: 

(13) 

For a given dual spin configuration {r^}, as shown in Fig. | we can draw 
a boundary one-cycle 70 on the original lattice, such that all the dual spins 
Tm inside 70 have the same sign. This is a two-to-one correspondence. In 
other words, for a given boundary one- cycle, one can assign either positive 
spins or negative spins inside it. For {r^} we get similarly a correspondence 
to 7q. Then we have 



2E E n 



-Jd 



I 



n 

7eCi ^^j^^j-^gp^ <mn>*G7+7o <mn>*^7+7o 

n e"^^ n n ^-^^ n ^""^ 

<mn>*e7+7o <mn>*^7+7,'| *697 



n 



n 



n 



n 



ymG<5(7o+7o) m^5(70+7o) »7iG<5(7o+7o) ™^5(7o+7o) 

After renaming 7 + 70 as 7 and 70 + 7o as 70, 



(14) 
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'^D=2^Y. E n e-^^ n 

7GCI 7Q6ro <mn>*e7 <mn>*^7 

n e-^^ n e^^ne-^^Re^^ 

<mr!,>* G7+70 <r?in>*^7+7o je97 1^97 I 

n e^^- n n n ^-'^ 

^mG(57o m^S-yo m&S^o mfS-yo 

This is the low-temperature representation of Zo- Comparing (|12D and (Jl 
we can easily arrive at the dual relations 



3 Some limiting cases 

For some special limiting cases the model reduces to known ones which are 
solvable. Then the dual relations in (^0) may be regarded as the dual rela- 
tions between two different models. Sometimes the reduction of the model 
is non-trivial and interesting. 

3.1 The limit A ^ 

From eq.(^ we have Pd — > oo. All the dual plaquettes will be frozen to 
Vi = 1. This is a constraint for the dual gauge fields {l^„}. In the language 
of cohomololgy theory, {Vi} are two-forms, gauge fields {Knn} are one-forms. 
Due to the constraint, {Kran} reduce to one-cocycles. 
The exact one-cocycles can be written as 

ViT'''' = VmVn (16) 

where r]m,f]m = ^ are also spin variables. That is, {Vi^r''^^ are pure gauges. 
It is known that all one-cocycles can be classified into a = 2^^ cohomology 
classes by the equivalent relation: and {Knn} are cohomologically 

equivalent if {Kn„V^„} is an exact one-cocycle; otherwise, they are not co- 
homologically equivalent. 

In our case, the cohomological equivalence is just the same as the local 
gauge equivalence. Cohomology classes are trajectories of gauge fields under 
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the local gauge transformation. If one can find an arbitrarily specified one- 
cocycle for each cohomology class, say {V^„}, then the general solution of 
the constraint is 

Vran = Vj^nVmVn, = 0, 1, . . . , a - 1 (17) 

The specified one-cocycle of the k-th cohomology class can be constructed 
from a one-cycle in the k-th homology class. Let 7^ be an arbitrarily specified 
one-cycle of the k-th homology class on the original lattice, one can define 

vL = v„M^\-' (18) 

mn [ 1 < mn > f. 7fc ^ ' 

Therefore the dual partition function in the limit A ^ is 



(19) 

Here the factor 1/2 arises due to the fact that {rjm} and {—rjm} correspond 
to the same gauge field Vmn- One can redefine the spin variables TmVm Tm 
and T^rjm — > t^, then the dual relation Z = c Zd, (||), reduces to 



^ exp ^ {m,s,s, + J'U,,S-S-) + i3uA = 

{Si,S[,U,j} [<ij>eL^ mGL^ J 

C II II exp i Vmnilk) {joTmTn + J'dV^ 

+ m AdT^T^ > 

(20) 

with c = 22^-i--P/2(sinh2Jsinh 2j')^/^{smh2p)P/'^ .The left-hand side is a 
'massless' two-spin system with gauge interaction. The right-hand side is its 
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dual model, whose partition function contains several terms due to the non- 
trivial topology of the lattice, of which each describes a 'massive' two-spin 
system with an anti-ferromagnetic chain. Both models are very interesting. 
The techniques of matrix representation may also be applied for solving them. 

For the planar lattice we have a = 1 {g = 0). When Ad — * oo, the dual 
model reduces to a simple Ising model and the phase transition appears at 
Jd + Jd ~ J* ■ This is a critical line. If one comes back the original model, 
the critical points lie on the curve 1/2 In coth J + 1/2 In coth J = J* when 
(3 also goes to zero. 



3.2 The limit A ^ 0, (5 ^ oo 

Let us further reduce the model in last subsection by taking the limit (3 oo, 
which implies A = 0, see eq.(|^). Now the model becomes self-dual and the 
self- dual relation is 



E E expj E u,,{^^){js,s, + j's[s'^i 

Q 1 I 

CE E 6^P) E ymn{lk)[jD 

where c = 2^"-^ (sinh2Jsinh 2j')^/^ 

and 7^ is an arbitrarily specified one-cycle in the /c-th homology class of the 
dual lattice. 

For the planar lattice the model breaks down into two decoupled simple 
Ising models. The critical points are two lines J = J* and J = J* in the 
{J, J ) plane. 




3.3 The Umit / ^ 

Here the spin 5*^ are decoupled. For the dual model j'^ ^ oo (or Jd — > oo). 
Therefore KnnTmT„ = 1, i.e. Vmn = 'Tm'''n- Now {Knn} are pure gauges and 
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oo 

i In coth J + i In coth J' = 



Figure 4: This is a phase diagram for a planar lattice. Here the real line 
corresponds to critical points at A = 0, the bold real line corresponds to that 
at A = cxD, and the squares are those for arbitrary A. The phase structure 
inside the box is unknown. 

Vi = 1. Therefore the dual relation (^) becomes 
^ exp < Y JUijSiSj + Y 

= C ^ exp <^ Y JoTmrn + ^Dr,n 

with c = 2^+^/2-^/2(sinh 2 J)^/2(sinh 2/3)^/2 rj^^jg jg known dual relation 
discussed by Wegner Q and Balian et al 0. 

In this case when /3 oo a phase transition occurs at J = J*. But note 
that now A is arbitrary. Correspondingly when Ad = and Jjj ^ oo there 
are critical points at Jd = J* for arbitrary [3d- 

Including two trivial cases, i.e. J = J =0, /5— s>oo and A = 0, J — oo, 
J ^ oo , all the critical points discussed above for a planar lattice are 
plotted in Fig ^. According to the analysis of the mean field method, inside 
the box there should be a critical surface. How to locate this critical surface 
analytically or numerically is at present under consideration. 

4 Discussion 

We close the paper by adding some remarks on the Wilson loop. 

In the limit A ^ and (3 —>■ oo, the Wilson loop along an one-cycle u) ( 
when cij is a non-boundary one-cycle the Wilson loop is also called Polyakov 
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71 



71 



Figure 5: This is a regular torus. 71 and 71 devide the torus into two parts, 
and inside each part the spins have the same sign but the two parts are in 
opposite signs. 



string) is 



< n u.^>- 



1 / ~ , ~ , ,^ (24) 

I] exp {JDVmn{lk)rmrn + JDVmn{lk + Uj)T^T^j 

^ fc=0 {rm,rL} <mn>&L\^ 

The Wilson loop is a topological invariant, i.e. its value only depends on the 
homological class of uo. This is because when a boundary one-cycle is added 
to cj, it can be absorbed into {t^} and the result is unchanged. Particularly, 
when is a boundary one-cycle, the Wilson loop is equal to one. 
Let us further take the limit J — » (J^ 00). Then 

< n >= ^ E exp E JoyLi^yJn (25) 

Consider a Li x L2 regular square lattice with periodic boundary conditions. 
The lattice is topologically equivalent to a torus, i.e. (7 = 1. 

In the limit — 0, everything is frozen and < n<ij>Gc^ > is trivially 
equal to one. However, suppose is very large, for a Polyakov string u; = 71 
the configuration shown in Fig ^ is dominating, where 71 can freely move. 
Therefore 

< n f/^, >~ ^le-'''^-'^ (26) 

<ij>e7l 

This implies the existence of topological excitations. 
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